On the Relativistic Foundations of Pseudospin Symmetry 

in Nuclei 



00 



o 
o 



J.N. Ginocchio^'* and A. Leviatan^'^'''' 
^ Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 
Racah Institute of Physics, The Hebrew University, Jerusalem 91904, Israel 



X) ■ Abstract 

We show that the generators of pseudospin symmetry are the non - relativis- 
ts^ I tic hmit of the generators of an SU(2) symmetry which leaves invariant the 
0^ . Dirac Hamiltonian with scalar and vector potentials equal in magnitude but 



opposite in sign, Vy = —Vs- Furthermore, within this framework, we demon- 



I strate that this symmetry may be approximately conserved for realistic scalar 

and vector potentials. 
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Pseudospin doublets were introduced almost thirty years ago into nuclear physics to 
accommodate an observed near degeneracy of certain normal - parity shell-model orbitals 
with non - relativistic quantum numbers (n^, i,j = i + 1/2) and {rir — l,i + 2,j = i + 3/2) 
where Ur, i, and j are the single-nucleon radial, orbital, and total angular momentum 
quantum numbers, respectively The doublet structure is expressed in terms of a 

"pseudo" orbital angular momentum £ = i + 1 and "pseudo" spin, s = 1/2. For example, 
{nrSi/2, {rir — l)c?3/2) will havc i = 1 , {nrPz/2, i^r — l)/5/2) will have i = 2, etc. These 
doublets are almost degenerate with respect to pseudospin, since j = £ ± s for the two 
states in the doublet. This symmetry has been used to explain features of deformed nuclei 
0], including superdeformation [^] and identical bands 0J^, and to establish an effective 
shell - model coupling scheme 0. Therefore there is an interest in understanding the origin 
of this "symmetry" The work of |^ has shown that a microscopic normal— i>pseudo 



transformation leads to an effective reduction of the spin-orbit splitting in the single-particle 
energy spectra derived from realistic (one boson exchange) internucleon potentials, in line 
with relativistic mean- field estimates The work of |]TD| has shown that quasi-degenerate 
pseudospin doublets in nuclei arise from the near equality in the magnitudes of an attractive 
scalar, Vs, and a repulsive vector, Vy, relativistic mean fields, Vs ~ —Vy, in which the 
nucleons move. A near equality in the magnitude of mean fields seems to be a universal 
feature of relativistic theories ranging from relativistic field theories with interacting nucleons 
and mesons |[ll[], to nucleons interacting via Skyrme-type interactions [|l^,|l^, to QCD sum 



rules [M- Recently realistic relativistic mean fields were shown to exhibit this approximate 



pseudospin symmetry in both the energy spectra and wave functions |T5|. In this paper we 
identify the algebra responsible for pseudospin symmetry as a non - relativistic limit of an 
SU(2) symmetry of a Dirac Hamiltonian with Vs = —Vy. 

The Dirac Hamiltonian, H, with an external scalar, V5, and vector, Vy, potentials is 
given by: 

= a • p + /?(m + ^5) + VV , (1) 
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where we have set ^ = c = 1 and a, f3 are the usual Dirac matrices |T^. The Dirac 
Hamiltonian is invariant under an SU(2) algebra for two limits: V5 = Vy and Vs = —Vy 
p!7Hl9[, (V5, Vy are spin scalars). The former limit has application to the spectrum of 
mesons for which the spin - orbit splitting is small. We shall show that the latter limit has 
relevance to understanding the origin of pseudospin symmetry in nuclei. 

The generators for the SU(2) algebra, S'j, which commute with the Dirac Hamiltonian, 
[if , ] = 0, for the case when V5 = — Vy are given by 



^ _ Q • P g^ Q • P (1 + /^) ^ . (1 -/?) f^s. 

where Si = ail2 are the usual spin generators and ai the Pauli matrices. This reduces to 
where 

^ 2 s • p 

Si = Up Si Up = — — Pi- Si . (4) 

In d^) Up = ^ is the momentum-helicity unitary operator introduced in [^] that ac- 
complishes the transformation from the normal shell model space to the pseudo shell model 
space while preserving rotational, parity, time - reversal, and translational invariance. What 
is significant is that the same unitary transformation appears in the non - relativistic limit 
of the generators of a symmetry possessed by the Dirac Hamiltonian for Vs = —Vy. In this 
case, the eigenf unctions of the Dirac Hamiltonian, H'^ = 8"$ are also doublets (5* = 1/2, 
/i = ±1/2) with respect to the SU(2) generators Si of Eq. (P|) 



S± vl>- = ^(1/2 T jl) (3/2 ± /2) v|/^±i , (Vs = -Vy) (5) 

where S± = Sx ± iSy. 

In general, the eigenfunctions of the Dirac Hamiltonian, have an upper and lower 
components, = -^^^^ ^I/. In nuclear spectroscopy, the lower component, which is 



small |T^, is usually ignored. The upper component satisfies the second order differential 
equation, 



(6) 



where 



H+ = HT + Vv + Vs 



(7) 



and the Hamiltonian is given by 



1 



£ + m + Vs-Vv 



(T • p . 



(8) 



Using the identity cr • per • p = p^, it is easy to show that 



Si , H. 



ps 



; 



(9) 



hence H^^ conserves pseudospin. As seen from Eq. the Sj are obtained by a unitary 
transformation from the ordinary-spin operators Sj, hence generate an SU(2) algebra 



Si , Sj 



(10) 



This means that, for Vs = —Vy, the non-relativistic single - nucleon wave functions (Dirac 
upper components) can then be labeled by pseudospin, s = 1/2, its projection, jl, jl = 
±1/2, and that these doublets are degenerate. The pseudospin generators connect 

wave functions of the doublets in the usual way. 



(l/2T/2)(3/2±/i) v[/+^^±i , 



'Vs = -Vv 



(11) 



where s± = s^ ±isy. 

However, in the exact pseudospin limit, there are no bound Dirac valence states for 
realistic mean fields [|I^] and therefore nuclei would not exist if pseudospin symmetry were 
exact. Nevertheless, it is possible to have small pseudospin symmetry breaking and have 



the requisite number of bound Dirac valence states 0,|T^] for nuclei to exist. Furthermore, 
we claim that the relationship between the wave functions of the doublets given in (|Tl]) will 
still be approximately valid. In order to prove this we see from (|^), that, in the pseudospin 
symmetry limit, the lower components are ordinary spin doublets with spin, s = 1/2, and 
projection fl, fl = ±1/2: 

A'' ^ —,fj, 1 



= ^(1/2 T Ai)(3/2 ± fi) ^_,^±i , {Vs = -Vv) . (12) 

These spin doublets are degenerate in energy and the wave functions "^-^p,, fl = ±1/2, have 



the same dependence on the spatial coordinates in the pseudospin limit. However, in [|T5 
it was shown that, even for realistic mean field Dirac Hamiltonians, the lower components 
of doublets have almost identical spatial wave functions, particularly for the doublets close 



to the Fermi sea. Even though |15| considered only spherically symmetric scalar and vector 
potentials, we assume this to be true in general and therefore label the states with spin pro- 
jection /2 referring to the dominant spin projection in the wave function. We now show that 
the pseudospin relations between the wave functions given in ( pTf ) for the upper components 
are approximately valid even when there is symmetry - breaking. 
From the Dirac equation, if^ = 8"$, we derive. 

Assuming that the lower component is normalized to unity, the normalization of the upper 
component is then (^+,^|^+,^) = Af'^ = (^_^^| q- • p )^ cr • p |^_^^). For the 

realistic Dirac equation with Vs 7^ —Vy, the lower component satisfies the second order 
differential equation. 



cr • p — ^— — a -p-m + Vv-Vs 

[t - m - Vs - Vv) 



|v[/_,^)=^|M/_,^) . (14) 



Using (|T^), we find that 



Af-' = {<iI^^f,\{S-Vv + Vs + m)^{S-Vv + Vs + m)\^_,f,) . (15) 



The normalized upper component will then be = Using similar manipulations, 

we then can derive the matrix elements of the pseudospin operators to be, 

($+,^,1 k > = Ar^{^-,ii'\{S -Vv + Vs + m)^{S -Vv + Vs + 

^ , (16) 

where the last step follows from the fact that the lower component has the spin and spatial 
parts approximately separated []T3| and from (|r^). Therefore, the near identity of the spatial 



part of the lower component wave function for the pseudospin doublet insures that the 
pseudospin doublets are approximately connected by the pseudospin generators even for 
realistic situations in which pseudospin is broken: 



1$+,^) ^ ^(1/2 T /i)(3/2 ± jl) 5^, , (Vs^-Vv). (17) 

In the case in which the potentials are spherically symmetric and satisfy Vs = — Vy, the 
Dirac Hamiltonian has an additional invariant SU(2) algebra; namely, 









(18) 



where £i = Up £i Up is the pseudo-orbital angular momentum operator, ii is the orbital angu- 
lar momentum operator, while ji = ii + Si = Up{ii + Si) Up = £i + Si. In this limit, the Dirac 
wave functions are eigenfunctions of the Casimir operator of this algebra, L • L 1$^^^,) = 
i{i + ^)\^ej m-), where we have used a coupled basis, j = i + s, where j is the eigenvalue of 
the total angular momentum operator Ji = Li + Si, J • J \^ijrnj) ~ + '^)\^ljmj)^ 
rrij is the eigenvalue of J^. Thus pseudo-orbital angular momentum as well as pseudospin 
are conserved in the spherical limit and 1^5 = — VV- (We note in passing that the operator 
K = —f3 {cr • i + 1^ which is conserved for arbitrary spherically symmetric Vs and VV po- 
tentials is given byi^ = 2 L • S + 1 ). From (p^ ) , we see that the lower component 
wave function will have spherical harmonic of rank I coupled to spin to give total angular 
momentum j. Since cr • p conserves the total angular momentum but p changes the orbital 



angular momentum by one unit because of parity conservation, Eq. ([T3|) tells us that the 
the upper component also has total angular momentum j, but orbital angular momentum 
i = i±l. If j = £ + 1/2, then it follows that £ = £ + 1, whereas if j = i - 1/2, then 
i = i — 1. This agrees with the pseudospin doublets originally observed |1|,0] and discussed 
at the beginning of this paper. However, the results here are very general and apply to non 
- spherical nuclei as well For example, for axially symmetric deformed nuclei, there is 

a U(l) generator corresponding to the pseudo-orbital angular momentum projection along 
the symmetry axis which is conserved in addition to the pseudospin for Vs = —Vy, 

where A = Up A Up. In this case the Dirac wave functions are eigenf unctions of A, A |$A,n) — 
M^Kn)' where Q is the total angular momentum projection, Q = A + fi, which has the same 
value for the upper and lower components since cr • p conserves the total angular momentum 
projection. Thus Q = A ± 1/2, corresponding exactly to the quantum numbers of the 
pseudospin doublets for axially deformed nuclei discussed in 0. Again we expect these 
multiplet relationships to be approximately valid for the realistic situation, ~ —Vy, 
following our previous discussion. 

The unitary transformation. Up (^, has been used to transform the non-relativistic 
spherical shell model space to the pseudo shell model space 0. Basically, as seen from Eq. 
(P!B|) and subsequent discussion, in the limit of = —Vy, this transformation transforms 
the normalized Dirac upper component wave functions into the normalized Dirac lower 
components. 

The Hamiltonian for the upper components in the realistic case which is given in 
leads to a Schrodinger equation. 



^ 'P^ + TT^ ^7 7r\^[(^ ■P,Vv~Vs](T -p + Vv + Vs 



(2m + Vs- Vv) (2m + Vs-Vx 



V 



= i£-m) 1$+,^) , (20) 
where we have made the approximation that (2m + Vs — Vy) » m — S, the binding 



energy. This Schrodinger equation is not of the form that follows from the usual shell 
model single - nucleon Hamiltonian because it has a spatial-dependent mass term and, since 
(Vv — Vs) /1m ^ 0.41 [O in the nuclear interior, the spatial dependence cannot be neglected. 



Furthermore, the second term in pO|), which produces a spin - orbit term for spherically sym- 
metric potentials, produces a momentum- dependent spin term as well. Although the helicity 
transformation of non-relativistic momentum-independent Hamiltonians and shell-model 



wave functions |21| induces momentum dependence, the above analysis demonstrates that 
the non-relativistic wave function (Dirac upper component) has momentum dependence be- 
fore the helicity transformation since it is a solution of a momentum-dependent Hamiltonian 
with a spatially dependent effective mass. 

One can ask if there are tests of the pseudospin doublet wave functions. One ingenious 
attempt was applied to polarized medium energy proton scattering for which the scalar 
optical potential is approximately equal and opposite in sign to the vector optical potential 

. Pseudospin predictions of the spin polarization and spin rotation functions were shown 
to be badly broken experimentally. The reason for this is as follows. In the strict symmetry 
limit the scattering functions for the partial waves depends only on £, not on the pseudospin 
and hence are equal for pseudospin doublets. When the symmetry is broken the difference in 
the scattering functions for doublets is proportional to 2£+ 1, just like the energy splitting for 
pseudospin doublets for bound orbitals ||T0| , p!5| . For medium energy proton scattering large 
orbital angular momentum, and hence, large pseudo - orbital angular momentum, dominate 
the scattering whereas in heavy nuclei the occupied orbitals have relatively small pseudo - 
orbital angular momenta (for example, the largest for an occupied bound orbital in ^''^Pb is 
i = 4). For this reason tests of the bound state wave functions may be more fruitful. 

In summary, we have shown that a relativistic symmetry of the Dirac Hamiltonian with 
Vy = —Vs reduces to the pseudospin symmetry in the non - relativistic limit. We have 
suggested that the pseudospin doublet wave functions of Dirac's upper components may 
still possess this symmetry even in the realistic Vy ~ —Vs situation. These results are 
expected to apply to non - spherical as well as spherical nuclei, in line with the validity 
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of the pseudospin concept shown for the oscillator shell model with arbitrary quadrupole 
deformations |2yj21|. We also conclude that the conventional non - relativistic shell model 
single-nucleon Hamiltonian could be made to be consistent with the relativistic mean field, 
and, as a consequence, pseudospin symmetry, by including a spatial - dependent mass and 
a momentum-dependent spin term in the mean field. 
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